arXiv:1508.04613vl [hep-ex] 19 Aug 2015 


PROCEEDINGS 

SCIENCE 



Measurement of Dijet Production in Diffractive 
Deep-lneiastic Scattering at HERA 


Stefan Schmitt*^ 

DESY, Notekestr. 85, 22607 Hamburg, Germany 
E-mail: 


sschmitt@mail.desy.de 


The production of dijets is measured in diffractive deep-inelastic scattering at HERA. The data 
were recorded with the HI detector at DESY in the years 2003-2007. Diffractive events are 
selected by requiring a gap in the rapidity distribution of the hadronic systen, where no particles 
are produced. Two jets are selected with transverse momenta in the hadronic-centre-of-mass 
system larger than 4 and 5.5 GeV, respectively. Cross sections are measured single- and double- 
differentially in various kinematic quantities. The data are found to be in good agreement with 
NLO QCD calculations based on diffractive parton densities determined from inclusive diffractive 
cross section measurements. 
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Figure 1: Diffractive dijet event as detected in the HI experiment. 


1. Introduction 

At HERA, reactions of electrons or positrons and protons, ep eX, are probed at centre- 
of-mass energies of 320GeV. At sufficiently large negative momentum transfer squared from the 
ingoing to the outgoing electron', > 4GeV, the process is referred to as deep-inelastic scattering 
(DIS), and defines a hard scale for perturbative QCD calculations. In fhe analysis presenfed here 
[||], dijef producfion in diffracfive DIS is sfudied. The reacfion may be written as ep —> eXY, where 
fhe sysfem X confains af leasf fwo hard jefs and fhe sysfem F is a proton or a low-mass excifafion, 
carrying a large momenfum fracfion of fhe incoming profon. A fypical evenf display is shown in 
figure ||. The low-mass sysfem Y escapes defection fo fhe leff (forward direction). The electron is 
detected in the right-most part of the detector. The system X consists of tracks and energy deposits 
in the calorimeter, forming two hard jets. There is a large gap in pseudorapidity^ between the 
acceptance limit of the detector in the forward direction and the most forward calorimeter deposit 
above noise level. 

QCD calculations for this process are based on diffractive parton densities (DPDFs). The 
DPDFs describe the probability to find a parfon wifh longifudinal momenfum fracfion zip in fhe 
profon, given fhaf fhere is a diffracfive signafure, characferised by a longifudinally momenfum 
fracfion 1 — x/p, momenfum fransfer t of fhe oufgoing profon and a hard scale p. They are folded 

'Throughout this paper, the term electron or the variable e is used to denote both electrons and positrons, unless 
otherwise stated. 

^The z axis is pointing along the proton flight direction. Polar angles 6 are measured with respect to the z axis. The 
pseudorapidity is defined as T] = — lntan(6/2). 
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with hard matrix elements to describe jet production at next-to-leading order in the strong coupling. 
Experimentally, the DPDFs are determined from inclusive diffractive measurements, where no 
requirements on the hadronic system X are made. As suggested by the scheme shown in figure 
the DPDFs are determined with the ad-hoc assumption that they factorise into a probability to find 
a colourless objecf IP in fhe profon and parfon density funcfions, ascribed fo fhe sfrucfure of IP. 
The probabilify or “flux factor” is faken fo depend on xip and t only, whereas fhe parfon densify 
funcfions of IP only depend on fhe variable zip and fhe hard scale /i. For fhe presenf analysis, fhe 
HI 2006 DPDF fil B 0 is used to predicf cross sections. Several measuremenfs af HERA have 
confirmed fhe validify of fhe NFO calculations in diffractive DIS i B1i- The presenf 
analysis is exploifing fhe high sfafisfics HERA II sample, in order fo obfain more precise resulfs 
and enable double-differential cross secfion measuremenfs. 


2. Data analysis 



The analysis |[T]] is based on evenfs recorded wifh fhe HI 
defector 0 in fhe years 2005-2007, corresponding fo an infe- 
grafed luminosity of 290pb^'. DIS evenfs are identified wifh 
an election in fhe HI rear calorimefer (SpaCal) [|g]. This 
limifs fhe accessible kinematic range in momenfum Iransfer to 
4 < <2'< lOOGeV^. The hadronic final sfafe X is reconstiucfed 
from fhe defected fracks and calorimefer deposifs, using an en¬ 
ergy flow algorifhm. The momenfum fransfer and fhe inelas- 
ficify y are reconstructed with the eX method [jll]], which uses 
properties of both the election and the hadronic final sfafe. The 
inelasficify is restiicted fo 0.1 <y <0.7. 

Diffractive evenfs are furfher selecfed wifh fhe condifion of 
a Farge Rapidity Gap (ERG) separating fhe hadronic systems X 
and Y. The ERG condifion requires fhat fhe most forward energy 
deposit above noise level in the HI calorimeter has a pseudora¬ 
pidity rjmax <3.2 and that there are no signals in the HI forward 
muon system or in the forward tagging system. The latter devices 

extend the reach for vetoing particles up to T] ~ 7.5. The momentum fraction xjp = {Q^ + All)/ i^y) 
is restiicted to the range xip < 0.03. Here, Mx is the invariant mass of the hadronic system X, 
reconstructed from all detected particles after excluding the election. The system Y is not de¬ 
tected; however from energy-momentum conservation one can deduce that its mass is limited by 
the diffractive selection to My ^1-6 GeV and the momentum transfer at the proton vertex is limited 
to |t| < IGeV^. 

The selection criteria are summarised in table The final sfafe X is boosted fo fhe y*p 


Figure 2: Diffractive dijet pro¬ 
duction in deep-inelastic scatter¬ 
ing at HERA. 


restframe and jets are reconstructed using the inclusive kp jet algorithm [12| with Pj recombination 
scheme and distance parameter /? = 1. Dijet events are accepted if there are at least two jets 
with transverse momenta of the leading (subleading) jet fulfilling fhe condifions pj j > 5.5GeV 
(P7-,2 > 4GeV). The jef pseudorapidify is restiicted in fhe laborafory frame fo —1 < < 2, in 

order fo ensure fhaf fhe jefs are well confained. 
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Extended Analysis Phase Space 

Measurement Cross Section Phase Space 

DIS 

3<Q^< lOOGeV^ 

4<Q^ < lOOGeV^ 


y <0.7 

0.1 <y <0.7 

Diffraction 

X[p < 0.04 

xjp < 0.03 


ERG requirements 

|f| < IGeV^ 

My < 1.6GeV 

Dijets 

pfi > 3GeV 

j >3GeV 


i >5.5GeV 

i >4GeV 


-2 << 2 '<2 

-1 < t]}"5 <2 


Table 1: Analysis phase space. 



The events are selected in an ex¬ 
tended analysis phase space, in order 
to ensure that migrations near the phase 
space boundaries can be well controlled. 

After correcting for detector effects, the 
phase space is restricted to the bound¬ 
aries described in the text. For unfold¬ 
ing from detector objects to the particle 
level, a regularised unfolding technique is 
applied. Typically, the distributions are 
measured in twice as many bins as are un¬ 
folded later. Choosing more bins recon¬ 
structed than measured later, improves the 
unfolding results and reduces the corre¬ 
lations between data points. The result¬ 
ing correlation coefficients are typically 
close to zero and do not exceed the range 
\pij\ < 0.6. In fact, only for a few cases, 

where the detector resolution is limited, correlations coefficients |p,yj close to 0.6 are observed. 
In addition to the fine resolufion on defecfor level, fhere are exfra bins ascribed fo confrol regions 
oufside fhe measuremenf phase space. The extra bins are also unfolded such that the normalisation 
of the prediction outside measurement phase space is taken from data. This procedure enhances the 
stability of the analysis to systematic effects related to the RAPGAP prediction and its reweighting. 

Control distributions of reconstructed variables at detector level are shown in figure |[ These 
are fhe momenfum fransfer (^, fhe leading jef fransverse momenfum, and fhe diffracfive momenfum 
fracfions xjp and zip, where zip = {Q^ / {Q^ +M\) and M 12 is fhe invarianf mass of fhe dijef 

sysfem. The variables are all well described by fhe prediction, which is used fo describe migrafion 
effecfs befween particle level and defecfor objecfs for fhe unfolding procedure. The prediction is 
evaluafed using fhe RAPGAP evenf generator and a GEANT-based simulation of fhe HI defecfor. 
If has been reweighfed in fhree variables fo improve fhe description; figure ^ shows fhe variables 
affer reweighting. 


Figure 3: Control distributions of selected kinematic 
variables at detector level. 
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Figure 4: Single-differential cross sections for diffractive dijet production in DIS as a function of the vari¬ 
ables Q^, y, xip and zip. 
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Figure 5: Single-differential cross sections for diffractive dijet production in DIS as a function of the vari¬ 
ables pj^, pj 2 , (pt) and At]*. 


3. Single-differential cross sections 


Single-differential cross sections for diffractive dijet production in DIS are measured in the 
phase-space defined in table ||. The dominant systematic uncertainties are: the hadronic energy 
scale uncertainty of 1%, model uncertainties of the RAPGAP prediction and normalisation uncer¬ 
tainties related to the LRG selection. The normalisation uncertainties amount to about 7%, the 
model uncertainties have a size of typically 5% and the hadronic energy scale causes uncertainties 
on the cross section of typically 4%. 

Single differential cross sections are measured as a function of the variables Q^, y, logxjp, 
Zip, Pt.I’ Pt,2, the average transverse jet momentum (pp) = {pj j + Pj2)l'^ and the difference in 
pseudorapidity At]* = |t]i —1]2|- The results are shown in figure ^ and ^ The data points are shown 
with their statistical and total uncertainties. The precision of the cross section measurements is 
limited by systematic effects in most cases. They are compared to NLO QCD predictions, obtained 
using the NLOJET-i-i- program [|^, 14] The predictions are shown with their DPDF uncertainties 
and their total uncertainties, then also including scale uncertainties. In general, the data precision is 
better than the precision of the prediction. In most cases the data precision is such that it is superior 
to the DPDF uncertainties, so the data have the potential to improve the DPDF fits. This is most 
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Figure 6: Double-differential cross sections for diffractive dijet production in DIS as a function of the 
variables and zip and as a function of the variables and pj ^. 


evident for the zip distribution. It is worth to note that no data with zip > 0.8 has been included in 
the determination of the DPDF fit at the time. However, because the NLO predictions suffer from 
sizable scale uncertainties, it is difficult to draw further conclusions or to pursue a DPDF fit, unless 
NNLO calculations for jet production in DIS become available. 

4. Double-differential cross sections and extraction of the strong coupling 

Double-differential cross sections are measured as a function of the variables and zip and 
as a function of the variables and pj j. The results are shown in figure ^ Similar to the case of 
single-differential cross sections, the data precision is not limited by the statistical errors in most 
cases. The data are well described by the NLO QCD predictions, with the exception of one point 
in the {Q^,Zip) distribution. At high zip > 0.5 the data are more precise than the prediction for 
any choice of Q^. In order to evaluate the sensitivity of the data to parameters of the QCD pre¬ 
diction, a fit of the strong coupling constant is performed, using the cross sections measured 
double-differentially in and pj j. The fit yields a^(Mz) = 0.119±0.004(cxp) ±0.002{had) ± 
0.005(DPDF) ±0.010(/rr) ±0.004(/r/)(ll). This first determination of in hard diffraction at 
HERA is caiTied out to as a consistency check and to evaluate the level of experimental as compared 
to theoretical uncertainties. The experimental precision of ±0.004 is better than the uncertainties 
stemming from the knowledge of DPDFs ±0.005. The largest uncertainties are related to limita¬ 
tions of the NLO calculation, ±0.010 and ±0.004 from the variation of the renormalisation and 
factorisation scales, respectively. 
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5. Summary 

Diffractive dijet production in DIS is measured at HERA, using data recorded with the HI 
detector. As compared to previous measurements on smaller data samples, the precision has been 
improved significantly, and is now limited by systematic effects. Single- and double-differential 
cross sections are determined for a selection of variables. The data are consistent with NLO QCD 
calculations based on DPDF fits from inclusive diffractive data. This confirms fhaf factorisation 
holds in diffraclive deep-inelasfic scaffering. Using fhe NLO calculafions, fhe sfrong coupling is 
exfracfed for fhe firsl time in hard diffraclion af HERA, and found to be consisfenf wifh ofher Us 
measuremenfs. The NLO QCD predictions, however, suffer from large scale uncertainties. Furfher 
inferprefafion of fhe dafa will benefil from NNLO calculations of Jef production in deep-inelasfic 
scaffering. 
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